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SUMMARY 


Investigation of the turhulent-houndary- layer flow over a flat 
plate in compressible flow is carried out on the basis of the scheme 
established in NACA TN 25^2. By averaging the Navier-Stokes equations, 
differential equations for the mean flow are obtained. A temperature- 
velocity relation follows without a specified form of the length scale. 
To derive the velocity distribution in the boun d ary layer, a choice of 
the length scale has to be made. The temperature-velocity relation 
reduces to Reynolds' analogy and the velocity distribution goes back to 
Von knrmnri's logarithmic law for the special case of incompressible 
flow. 


There are essentially three universal constants, arising out of the 
correlations in the energy equation, to be determined by comparing with 
suitable experiments of the teii 5 >erature- velocity relation at any known 
Mach number and heat transfer at wall. The behavior at other Mach 
numbers and heat-transfer conditions may then be readily predicted. 
Because of the lack of acctirate experimental data, attempts to carry 
out such determinations are not included in the present report. 


INTRODUCTION 


As part of an investigation of Von Karmen' s similarity theoiy and 
its extension to compressible flows, the theory for incompressible flows 
was examined critically in reference 1 by tising modem concepts. It was 
found that the original form of the theory is supported by modem con- 
cepts. In reference 2, the theory was extended to the case of compres- 
sible flows and it was found that the analysis could be carried through 
with additional approximations but without any modification of the basic 
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concepts. In the present paper the problem of compressible flow in a 
turbulent boundary layer is considered on the basis of the scheme 
established in reference 2. It is shown how the influence of the Mach 
nijmber (e.g. , on the velocity and temperature distribution) can be pre- 
dicted from the theory after the constant coefficients in'^the theory 
are determined by one set of experimental measurements. 

This investigation was conducted at the Massachvisetts Institute of 
Technology under the sponsorship and with the financial assistance of 
the National Advisory Committee for Aeronautics. 


VELOCITY- TEMPERATURE RELATION 


Attempts will now be made to deduce relationships among the meein 
quantities and their distributions within the boundary layer. To do so 
one returns back to the complete equations of motion (equations (26) 
to (29) of reference 2) instead of working with the "localized" equa- 
tions where the observer rides with the mean velocity. These equations 
are averaged so that only mean qiiantities and the correlations of the 
fluctuating quantities appear. 

A strictly parallel flow will be considered. The Justification of 
its application to the case of a boundary layer will be discuased later. 
By using the continuity equation, it is easy to show that 



( 1 ) 


when the mean flow depends on y only. Here the bar denotes the average 
of the quantity under it and F stands for any function. (See appendix A 
for definitions of symbols.) Thus, equation (26) of reference 2 leads to 


^ puv = 0 





^y. 


— pvw = 0 

^y 


(2) 
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There follows immediately, by integration. 


puv = Constant = Tq 



( 3 ) 


pvw = Constant = ' 


where Tq is the shearing stress at the plate in the x-direction, 

Tq2) "tlie shearing stress at the plate in the z-direction, and p^, the 

pressure outside the boundary layer- One expands the averaged quan- 
tlMes in equations ( 3 ) and finds that triple correlations are involved 
as well as the mean-flow variables and double correlations. In usual 
cases the triple correlations are taken to be of higher order than the 
double correlations. For example. 


p'u' V* 


o|^(p' u* v' I 





u’v’ )/ « 


p u’v’ 


( 4 ) 


A word of caution, however, may be in order here to warn against the 
overconfidence in such estimations. It is nothing more than a plausible 
guess of the individual terms, and sometimes a conibinatlon of the indl- 
vidiial terms may very well Invalidate the resxxlt. The discussion in 
reference 2 of the dilatation e serves as a good example. 

Nevertheless, it is perhaps admissible in the present case to drop 
triple correlatlor^ involving p' on the basis of the estimations like 
equation ( 4 ) . If one concedes that the effect, of the triple correla- 
tions is not negligible but relatively small, then the validity of the 
theory, being an approximate one, might not be serioxisly impaired 
because of the omission. By leaving out the triple correlations, equa- 
tions (3) are expanded into the following: 
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p U*v' + U p’v’ = Tq 


(5) 


p = Pi + p (v* )2 




Pi 



( 6 ) 


p V*W’ = Tq2 (7) 

where subscript 1 denotes the free-stream quantity. From the equation 
of continuity, one has 


hence 



0 


p’v* = Constant = m 


( 8 ) 


where m is evidently the mass transfer through the wall. Likewise from 
the equation of energy (equation (28) of reference 2) 






which yields, by omitting the triple correlation p’v'T' and inte- 
grating between 0 and y. 



v*T* 


+ T p'v’ 


j - R(i 



r- 

|ie dy + q + q* 


J6* 


(9) 


where p^ is the pressure at the wall (cf. equation (6)), 6*, the 

effective thickness of the laminar sublayer, q, the heat transfer at 
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'^ 6 * 

the wall, and q* = [le dy, the heat generated within the laminar 

Jo 

sublayer (cf. appendix B for discussion of q*). 

One may now make the following ohseDrvations: In equation (8), the 

turhiolent mass transfer m vanishes if there is no addition or sub- 
traction of mass through the plate. For a nonporous flat plate, then, 
equation (5) leads to 


p u'v’ = Tq ( 10 ) 

One may note that equation (lO) is also the basic equation in refer- 
ences 3 and h, where the role of p'v’ was not mentioned. In refer- 
ence 5; the tenn p'v* is kept and given no physical interpretation. 

With the previous assimption of small turbulence level 

(v* )^/u^ « 1, when the free- stream Mach number is of order imity, 
eqiiatlon (6) reduces to 

P sj p^ (11) 

as usually accepted for botindary- layer flows, and the term p - Pq drops 
out in equation (9). It is seen, however, that eqiiation (ll) ceases to 
be true if 


Ml 


X 


(v*)' 


0 ( 1 ) 


U': 


Equation (7) is merely the statement of constancy of the transverse 
shear in the z-direction. It does not influence the mean motion in the 
xy-plane. 

The next step is to introduce the similarity theory and represent 
the fluctuations in terms of the scales. One obtains from equation (lO) 



( 12 ) 
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by absorbing bbe correlation constant u.'v' into Vq, tbe similarity 

scale for velocity fluctuations. Similarly :^om equation (9)^ witb T' 
broken into and Tg' according to reference 2, 




Cpplv'Ti' + v'T2' 



_ 2 

“■iP^o ^ 
to 


+ q + q* 


(13) 


where aq is a correlation constant for the dissipation and use has 
been made of the relation 



pe 



P 


ijte/p avo^/to 

Introducing sceiles 0q ajid ©2 1"°^ Tq’ and Tg* ^ respectively, and 
absorbing the correlations into proportional constants, 02 


Cpp(v'Tq’ + v‘T 2') = PVo(o2lo ^ + a3Vo2^ 

Using equation (12) and the definition of t^ in equations (32) of 
reference 2, 


— 2 dy 

i8* U 


py 


da 


5* 


= aqTo(u - u*) 


with u* representing the velocity at the edge of the laminar sublayer. 
Hence equation (13) becomes 
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or 


^ — + a3"'^o “ - u*) + (l4) 

where again the relations (12) and equations (32) of reference 2 have 
heen used. 

The quantity Vq may he put in terms of T: 



(15) 


with 3 = —• Also, since q* = t^^u* 

introduce an "effective" velocity 


(see appendix B) , one may 


Ue = u - 



(16) 


Then the final form of the differential equation between T and u is 
reduced from equation (l4) into 


02 


^ + 033^?= o^Ue + q 


(17) 


Equation (17) is derived without a specific assumption on the length 
scale Zq. The constant o^ represents the heating due to dissipation, 

as a part of the external work t^uj 02 represents both the mixing 
phenomenon and the compressibility; and 03 reflects the combined effect 
of the compressibility and the dissipation. On a closer look at expres- 
sion (16), it is revealed that the factor (1 - — ) arising from the 

\ ^ 1 / 

heating effect is originated from the idealized situation of an entirely 
turbulent outer boundary layer with similarity and an inner laminar 
viscous sublayer. In the sublayer all the energy supplied by the 
external work must go into heat through dissipation, hence the constant 
unity appears. In the turbulent layer, part of the energy is trans- 
ferred or tiirned into turb\ilent energy, therefore only a portion 

is dissipated. If a transition region between the two idealized regions 
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had heen ass-umed, there would he no discontinuity in the picture hut 
a term of the order ^1 - would still he present in the heating 

because of dissipation. ^ 

Equation (l?) inay "be nondimensionalized hy writing 



where subscript 1 again refers to the free- stream value and is 

the temperature at the wall. Then there follows 



(AiU + 



( 19 ) 


^o see this, one could assume a linear transition of the dissi- 
pation parameter such that 


— _ ^ du 


qe = Tq 0 < y < ci6*, c^ < 1 


1 + 


u - Uc 

(ai - l)= ^ 


Uc2 “ 


To C^s* ^ y ^ Cgh*, Cg > 1 


= ^ 1^0 ^ y 

where Cj^B* ^ y CgB* is the transition region and u^ and u^ are 
the respective velocities at the ends of the transition. Then the term 


corresponding to ^1 - is ^1 - — ^"""2 — ~J 
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where 


Ai = 7R — 

-L 02 




^ ^ Sw' 


( 20 ) 


Equation (19) is the fundamental equation of the present theory for the 
variation of mean temperature with mean velocity. In the spirit of the 
local similarity, since the tiirbulence pattern is obtained by letting 
the observer ride with the mean local velocity, it seems quite clear 
that the Mach number effect at most would be an indirect one. The 
constants a^, and being essentially correlation constants, 

OTight to be nearly independent of the Mach number. The effect of Mach 
number on the mean temperatTxre-velocity relationship is largely in the 
first term on the right-hand side of equation (19) ^ where M^^ stands 

before the parenthesis. One may also expect the constants A 2 and Q 
to vary in some way with the Mach number throu^ the frictional coef- 
ficient and the heat-transfer coefficient c^.- However^ the varia- 

tions of Cf and Cq with Mach number are known to be rather slow. 

For an approximate theory, it might be sufficient to regard A^^ 

and Q adso as \miversal constants independent of the free-stream Mach 
number. Likewise, the constants A^ and Aq might further be taken 

to be Independent of the heat-transfer situation at the wall. 

There are consequently three parameters 02 ^ and to be 

introduced in the problem besides the boundary conditions of the mean 
flow at the wall and in the free stream. These parameters are related 
to the turb^llence mechanism emd must be empirically determined by com- 
paring with experimental data. Once determined from suitable experi- 
ments at low speeds with known heat transfer at the wall, they enable 
one to predict the behavior at any other Mach number or heat- transfer 
condition, provided, of course, that the turbulence level is still 
small «nd that the Mach nimiber is not excessive. 

In the following, equation (19) will be Integrated first for the 
case of subsonic flow and then for the general case. For the subsonic 
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case the solution can he expressed in an ascending series of M]_ , and 
a few terms would he sufficient to compare with low-speed tests to 
deduce the values of the ■universal constants. The integral for the 
genereil case is ■useful mainly for the prediction of supersonic ho^undary 
layer after the uni versal constants are de'bermlned. 

Case of suhsonic flow .- Aesvme the solution to he expanded in 
ascending powers of 


0 = 0^Qj = %^®(l) %^®(2) + • • • + + • • • (21) 


hence. 


p- = e(o)l/2 + + 






( 22 ) 


2n 


Substituting into equation (l9) and equating powers of , one obtains 

a set of differential equations for the functions ®(o)^ ®(l)^ ’ * 

0^^^. Thus, for 


M 


0. 

1 • 




Ml 


M- 


6. 


f!(2l- 0 

dh 




■18(2) _ 1 -l/2„ ^ 

-^-2V(0) ' ®(1) 


"^^(3) - 1 a n -l/2[s i fl -2 

- 2 V(0) 1^(2) ^ ®(0) ®(1) J 


(23) 
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and so forth. For boundary conditions one has for the incompressible 
case with no heat transfer ® = ®(o) “ ^ when li = U]_^ where 



This result holds for any Mach number; therefore. 


0(l) = 0 ^ 2 ) = • • . = 0 when u = U]_ 

With these conditions the integrals of equation (23) are found to be 


0(0) - 1 + 


Q(S - Ui) 


'( 1 ) 


= I Ajlf - tii^) t 0,3/2 . 


6(2) = 


2Q 


In 5/2 , 

\ 2Ai(l - Q) 1 

3/2 A 

■ " 

) 3q2 ( 

,"(0) - "■] 


> (24) 


Aid - Q)^ 

L 


- Ai 


kkr 


- 1) * t|(«( 0 )^ - 1) 


and so forth. When the wall is insulated, Q = 0, and the functions 
0(l), ®(2).» • • • slnqjly pol 3 nQominals of u. For convenient 

reference, a few of them are listed as follows: 


«(o) - 1 

®(1) = ' ^ 1 ) 

®(2) “ - (t ^ *2)('a - 

•(3) - ¥ [■ - V) • (s - • 

[? - - '*‘1 
^ * ~2 ~ •*• 2 ^ 2 ^^ {% - + 

" f - " A3.A2 . AgS)] (S _ ^,)| 


and so forth. 
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An interesting result to be noted is that the solution ®(o) ^ 

equation (24) evidently coincides with what is usually known as Reynolds' 
analogy, where a similar transfer mechanism is assTimed for both the 
momentimi and the heat transfers. Here the same results are obtained 
because of the fact that the turbulent exchange, proportional to the 
temperature gradient dl/dy, dominates the situation at very small Mach 
numbers. A similar tr ans fer mechanism follows the assumed similar tur- 
bulence pattern. The solutions ®(i)# forth now give the cor- 

rection to Reynolds' analogy at higher Mach numbers. 

General case .- A brief outline will now be given for the general 
case. Rewrite equation (19): 


JdW AgMj^^X A^M^g ^ ^ g 

W 2 / 2 2 


(26) 


which is now homogeneous in and u. Using the standard method, put 


^ = F(u) 


[A^Mi^ 


u + l) 


(27) 


Then F is to be solved from 


F dF 


du 


W W 1 AiMi^ tg , Q 

2 2 2 2 


(28) 


Integrating equation (28), one gets 


+ I = c(f - ki)^l(F - kg)' 


(29) 
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where 


- 

= I 

Si = , 

®2 = ^/(^i - ^) 

C = Integration constant 



(30) 


Together, equations (2?) and (29) form a parametric representation of 
the function 6 = 0(tt). Applying the condition at the outer edge, 
namely, 0 = 1 when u = uj, one has 


and 



The complicated way in which the constants are entangled renders their 
evaluation from en5)lrical data highly tedious. Besides, the proper form 
of equation (29), where every term is real, depends on the sign and 
magnitude of the constants and is hest to he directly integrated from 
equation (28) once the constants are known. 


VELOCITY AND TEMPERATURE DISTRTBDTIONS IN 
BOUNDARY LAYER 


To determine the distribution of the mean velocity and/or the mean 
temperatTore across the houndary-layer thickness, it is necessary to bring 
in the definition of the length sdale Zq as a function of y. In 

reference 2, two alternative expressions (equations (33) and (39)) have 
been given for the length scale, laying eu^Jhasis on the momentum and the 
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energy equations, respectively. Integrating either vill lead to a 
velocity distribution, from -wMch a temperature distribution may be 
calciilated by means of the results in the section "Velocity- Temperature 
Relation. " The distributions evaluated from the two expressions irill, 
in general, be different. It seems that only by comparison with 
experiments can one say whether both lead to essentially the same 
result or one form is preferable to the other in certain particular 
cases. 

The two definitions are, nevertheless, of the same general form: 


^o 


^ dx/d% 
dy/dy2 


( 32 ) 


with X rtanding for the variable u in equation (33) of reference 2 
and for T in equation (39) of reference 2, Without specifying X, 
it follows from equation ( 32 ) of reference 2 that 


^o 


du 

dy 


= an 


du dX /d% 


1 dy dy/^y2 

where a^ is a proportional constant. Rewriting, 


(33) 


dX\ du 


Hence, by integration and substitution of equation (l5). 


log. 


dy 


du 


fl / ^ 

P J Tt 


( 34 ) 
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The temperatiore-velocity relation may then he used, yielding 


S » G(S) 

<iy 


where 



(35) 


Finally, one may obtain the distribution by integrating equation (35) 



'dX de du 


J d0 du G(u) 

For definition (33) of reference 2, X = u. 


y = 


du 

G(u) 


for definition (39) of reference 2, X = 0, 


(36) 


(37) 


y = 


n 


u 


du G(u) 


(38) 


Again one may first solve the subsonic case for a more explicit expres- 
sion, which mxist be reducible to the usual logarithmic distribution for 
incompressible flow. Then a brief discussion of the general case for any 
Mach niimber will be taken up. 


Case of subsonic flow .- The fxmctlon G(u) must first be evaltiated 
before integrating equation (36). By definition (35); 


G(u) 




= exp 


(39) 
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after nondimensionaliziiig with equation (8). Using equation (22), 
one has 


where 


'u 


^ = (Pq(u) + M^^cpi(S) + + 

u^ 


CPq(^) = 


Ou 


^1 




(Pl(u) = - 2 


■u 


-3/2 


(0) ‘^(l) 


6/tN (iu 




V m 


92^^) = 


'u 


Un 


(§ S(1)^«(0) 


-5/2 




du 


93 ( 1 ^) = 




li-. 


fl O fl -5/2 1 g, -3/2 


r '^(i)'^(2)"(o) 


-1 <iu 


I 6 "(1) (0) 


') 


and so forth, with the defined by equation (24) with heat 

transfer at wall «nd by equation (25) with an insulated wall. Sub- 
stituting into equation (39) ^ one may therefore put 
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G(u) = 


go(u) gi(u) + Mj^^ggCu) + . . .J 


n 


where 


go(u) = exp (Pq(u) + 




y (^ 1 ) 


g2(vi) = 


c?2(^) + ^ 


2 2 


(P-f (u) 


= ^1 \|^ 


<P3(u) a^\j 


4 <Pi(u)92(S) + ^ 1>i3(B) 


Inverting^ 


G(u) go(^)L ^ -1 ^ -2 J 


where 


(te) 


S_2^(n) = -g^(u) 


g_2(u) = g^^(u) - gg(u) 
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and so forth. Equation (42) may now he used to evaluate eqiiatlons (37) 
and ( 38 ): 

( 1 ) Case (a) X = u: 


pd 


V - Yr = 




du 


where y^. is a reference station, u^ is the value of u 
at y^, and 


u,0 


'u,l = 


^.2 = 


'u 


du 


u, 


gl(g) 


du 


Pu 


Ur 


gl^(ii) - g2(i^) 


du 


(43) 


It will he of interest to apply equation (4-3) to the incompressihle 
case, namely 


2 — ^ 0, The distribution is then 


y - y = Y „(u) 

^ U,0'' 
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For simplicity, one may take yj. = 5, 6 being the thickness of the 

p 

boundary layer, so that 





II _ til = ^ _ 
^1 


1 


With given by equation (24), 



^Actually this procedure does not give the best fit to ej^erlmental 
velocity distributions, since the similarity concept is not valid near 
the outer edge of the boxindary layer. Instead of 5, the length scale 
should be taken as proportional to v vhere u-^ is the frictional 
velocity. For general discussions, however, there ought to be little 
difference. 
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vhere 



The constant ao may be determined from the usual condition at the 

3 

vail. 


du 

dy 


00 


as 


0 


The other constant a^ has to he evaluated hy matching vith experi- 
mental distribution. From the definition (33)^ it corresponds to 
Von Karma’s universal constant. 


One may further develop the solution (44) into a power series in 
Q for cases where the heat transfer is relatively small. After some 
manipulation, the final result is 


where 


loge 




= -BiX * x) 




( 1 ^ 5 ) 



B2 = 





^As a possible refinement, the laminar sublayer may be assumed to 
have a linear velocity profile, and, instead of infinity, the velocity 
gradient at the wall may be prescribed as 



y 


0 


However, this step complicates the practical calculation and must be 
justified by experience. 
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Applying the condition of infinite velocity gradient at the wall, one 
finds easily 


Bg = -5 

For the case of an insulated wall, there follows 


(46) 


u 



(47) 


which is precisely a form of the logarithmic law, matching with the 
experimental distrihution at the assumed edge of the boundary layer. 
Such a form has been suggested by Dryden (reference 6) in 1935^ where 
the frictional velocity and Von Karma’s constant K were used 

the definition 


(48) 


Eq\iation (47) then reduces exactly to Dryden's form 

u-uJ = -^logg| (49) 

One may also conclude that the logarithmic law, suitably modified, 
indeed can be approximately true even with the presence of a very small 
amount of heat transfer, since equation (45) shows that the first-order 
correction for Q = 0 merely amounts to replacing the factor in equa- 
tion (49) by a slightly different one. If the heat transfer is appreci- 
able, equation (45) predicts that the logarithmic law in the form of 
equation (49) would break down. There would be considerable interest 
if equation (45) could be tested for different values of Q. 

It is apparent that each of the fractions T in equation (43) 

may be developed into a power series in Q as above. In general, 
therefore, the right-hand side of equation (46), which gives the velocity 
distribution, is in the form of a double series in the parameters Q 

and involving constants a^^ and 82 which are to be determined 

once and for all. In fact, by comparing with the incon5>resslble 


instead of a^ and c^ in equation (47). 
one verifies readily that 


“-“'ll 

ai = K 
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ins\ilated case as atove, "both a^^ and a 2 are given in terms of 
and '£> 2 - 

For its practical importance, the results for the inSTilated case 
at subsonic Mach numbers are ejcplicitly given as follows: Equations (4o) 

are integrated by using equations ( 25 ), 


<Pq(u) = u - Uj_ 

cpj^(u) = bQ + bju + + ^3^^ 

^ 


( 50 ) 


and so forth. 

where 

-1 Ap ^ 3 

A2 ~ ; 

■bo - - -^ ^1 

- T^i 

Aj ^ g 

Ag ^ 

hi = -jj- u^ + 

T"i 

I 

II 

<M 


At 

b, A 


3 12 


c^--— a 2 u 5 _J_ 

=0 - 20 H "1 16 


<^1 = 


+ h ^1^2^!^ + \ A^Ag + i 

~ ^ A^AgUi^ - I ^ A^Ag - I Ag^ 


J. 

32 


( 51 ) 


eg = 

^=3 = - ^ Ai'^ui^ ^ ^ A^AgUi + 


2~ 2 1 


8 


12 


cl^ = ^ AiAg 

c . 3 2 

5 160 ^ 
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and so forth. The functions of equations (4l) hecome 

gQ(u) = exp [bi(u - Uj) + a^ 
g^(u) = B3^qp]^(u) 

g2(u) = + I B^^qp^^^Cu^ 

and so forth. Hence ^ 


^u,0 = 


{i - [:®i(^ - ^i)]j 

Yu,l = ®2 - I^l)j(dl0 + 

^<3-xo + '^11^1 ■'■ ■*■ ^13^1^)^ 

^-exp [j-Bi(2 - %) - (<^0 + <^211^ + ^ 2 ^^ + 

d22^^ ■*■ dgijU^ + d^^u^ + d2gu^^ + ^*^0 "*" *^21^1 

d22^i^ + + ^24%^ ■'■ ^25^1^ ■*■ ^6^1^)] 

and so forth, where 


Yu, 2 = ®2 


J 


d_ = h 


10 ‘'0 B 


hi 2bp 6h:j 
+ — + — — + — - 


1 Bn^ Bn 


2bp 6b q 

dll = bx + :^ + :^ 


Bi Bi2 


3bq 

dx2 - ^2 + 3^ 


^13 “ ^3 


(52) 


(53) 


(54) 
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(m. + n)I ^ 


The constants h™ are defined as follows: 


B-i 0 

ho = -cq + ho 


hi = -c^ + B3_hQh^ 


h2 ^ -C2 + ^(^1^ + SbQhg) 

hj = .C3 + Bi(i=i 1.2 + Dobj) ■ ( 56 ) 

ilj. = -oij + 

h^ = -C5 + Bih2h3 

>1 - -h 2 

hg - hs J 

With these evalxiations, the distribution for the insulated suhsonic case 
becomes finally. 


logg g - 1^ + Mi2fi(u) + M/f2(u) 


+ • • • 


fl(u) . - ^ d4u= - %■>) 


[S) = ^ (^2n - I XI - “l”) 

n=0 ' Z+m=n / 
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and so forth where B 2 is now determined again by the condition of 
infinite velocity gradient at wall, „ 


_S_ 

B 2 


= -1 + Ml 


YL 

n=0 


<ilnV 



Yi 


n=0 


+ . . . 


The other alternative for X will now be coMidered: 
(2) Case (b) X = 0: 

Instead of equation (43), the equation becomes 


Pu 


y - jj. = 




du G(u) 


Since, from equation (43), 


G(S) duZ_ 


n=0 


~ Y 
1 dS 


one may rewrite eqiiation ( 59 ) as 


( 58 ) 


(59) 


where 




^0,n 


pu 


UUr 


y- ^^U) *^yu,m 

zfen 


du 


y ( 60 ) 


In parallel with case (a), where X = u, the results for incompressible 
flow with heat transfer and for an Insiaated wall in subsonic flow are 
given in the following. For the former problem, by using equation (23) 
one obtains Immediately 


''0,0 = ^u ,0 


( 61 ) 



26 


MCA TN 2543 


60 that a modification of the definition of B2 leads to the same 
distribution (^5). Thus the two choices of X are equivalent to each 
other even for the distribution in incompressible flow. 


For the compressible case with an insulated wall, the right-hand 
side of equation (60) starts with Now 


d0 




(1) <lYu,0 


du 




du 


du 


u 


Ui 


(a^u + A2) 


dX 


u,0 

du 


du 


by taking y^ = 5 and using equations (23) and ( 25 ). Integrating, 


^e,i = 


^ + b7 ^1^) 


^-Bi(u-ui) 




Substituting into expression (59)^ one has 


y - 5 = + o(mi^) 

The condition of infinite velocity gradient at y = 0 again 
determines B2, 

6 = + h + AiuJ + o { m /) 

Hence the distribution may be written as 



(62) 
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Equation (62) indicates ttiat in the absence of heat transfer at the 
wall, the choice of X = 0 leads to an additional term to the usual 

logarithmic law even when « 1 (hut not zero). The additional 

term is of more weight nearer the wall where u^u^ differs appreciably 
from unity. Noticing further that A^jA^ is a controlling parameter, 

A2 

by equation (20), one may say that the predicted deviation from the 
visual logarithmic law depends on the relation between the heat gen- 
erated by dissipation and the turbulent transfer of such heat (and 
due to compressihllity) . Ordinary theories assume the dissipation to 
he negligible, equivalent to putting = 0. Then the additional 

term vanishes. 

The apparent contradiction of equations (45) and (62) when only 
the lowest -order terms in both are retained is resolved by the obser- 
vation that, in complete form, one should have 


y - 


5 = or, 


u,0 




( 63 ) 


so that the lowest-order terms should he taJcen according to whether 
Q « or the converse is true. If both are of the same order, 

neither should he left out and the final expression certainly -viill 
contain an additional term analogous to that in equation (62). 
Physically, it is obvious that if a great amount of heat conduction 
is present at the wall, it certainly would overshadow the distributed 
heating due to dissipation (and compressihllity) at lower Mach numbers. 
On the other hand, when there is no external heating, the dissipation 
becomes the only predominant factor. 

The Mach number effect in ecjuatlon (62) can he worked out without 
difficulty, but, involving no new features, the details are omitted in 
this report. 

General case .- The general case will he only briefly treated, as 
the integration depends on the values of the constants A-^, A£, and Q. 
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(Cf. equations (28) to (30).) Still starting with equations (36) 
and (39)^ one substitutes ^ hy equation (27) 


where 


G(u) = exp 


n~ 




du 




^1 F(u)U!^^ 




a2 


= exp -B-| 


dF 




_ 






F - 1 


= B 


h - ki\Pl 

V - '^1 


Pi = 




B^ = e 


a2 


f(ui) - k2" 


FI 


(“1) - 


Pi 


(64) 


Hence, for X = u. 


y - 


u 


du 



G(u) 

1 

Op 

du| 

®3 
J u 

F^dF 


(f - kiVPl 


dF 


(65) 
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■with du/dF given by equations (28) and (29). Similarly^ for X = 0, 


y - Yr 


®3 



d0 du p - ki VPl 
dxi dP - kg/ 


dP 


( 66 ) 


with d0/du to be evaluated from -equations (27) and (29). The con- 
stant of integration is determined from the condition: 



as 


0 


The constant has previously been identified with 

CORRELATION OP THEORY WITH EXPERIMENTS AND DISCUSSIONS 
Determination of Arbitrary Constants 


There are two sets of arbitrary constants. In the tenrperature- 
velocity relation, one has three constants cxg, and representing 

essentially combinations of the correlations, which are assinned uni- 
versal with regard to the free-stream Mach number and the heat-conduction 
conditions at wall. In addition there is the mean velocity u-*- at the , 
edge of the laminar sublayer sepeirating the turbulent boundary layer 
from the wall. The second set arises out of the Integration for the 
mean velocity distribution within the boundary layer, consisting of 
B^ and Bg in the section "Velocity and Temperature Distributions in 

Boundary Layer. " The external conditions defining the mean flow are 
the foUo'wlng: The free-stream velocity, pressure, temperature, the 

friction at wall, and the temperattire at wall (or the amount of heat 
conduction from the wall). 

With given external conditions, the constants can be determined 
by comparing •with accurate experimental results, G?he significant point 
is that, once determined, a \miversal law for all Mach numbers and heat 
conduction is ' established. No other ad hoc assumptions need to be 
introduced for special ceises. 
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It is then only necessary to conduct a few low-speed experiments 
with prescribed heat conduction at the wall and measure ^the mean tem- 
perature and velocities at points across the section. Eli^' data 
(reference 7) were in this category but unfortunately were not pre- 
sented in enough detail to suit the present purpose. One may nattirally 
also deduce from high-speed meastirements with somewhat more work. 

Wilson (reference 8) and Ladenburg and Bershader (reference 9) published 
the res-ults for a supersonic free stream. Unfortunately^ again, their 
data were not sufficient, because only one of the two variables u 
and T was measiared in either case. Both assumed the isoenergetic 
relation to hold for evalTiation of the other variable. 

There are arguments for the validity of the isoenergetic relation, 
either from the practical reason that the effective Prandtl number is 
nearly unity or from some assunrption on the mixing lengths for momentum 
and energy, such as Ferrari' s (reference 5). Even admitting the argu- 
ments, no isoenergetic law would follow when the dissipation term is 
kept in the energy equation. In reference 2, it has been shown that 
the dissipation terms should be retained for a consistent theory, and 
one could only regard the success, if any, to be essentially an 
empirical one for the Mach number or heat conduction involved. Such 
empirical results might be useful, however, for the approximate evalua- 
tion of the constants a]_, a^, and a^, at least as a guidance to the 

orders of magnitude. 

The constant u* is a more troublesome one. In appendix B, it 
is shovm to be dependent on the Mach number as well as another empirical 
constant defining the extent of the turbulent layer. The correct determi- 
nation, therefore, is very tedious and perhaps not warranted because 
of the approximate nature of the theory. A possible way is to deter- 
mine its veilue from incompressible flow and neglect the dependence on 
Mach number. The error Introduced can be assessed only by comparing 
with experimental, results. 

The determination of B]_ and B 2 is done by matching with the 

empirical distribution of the mean velocity or temperature. In the 
section "Velocity and Temperat'ure Distributions in Boundary Layer" it 
is found that B 2 is essentially determined by the condition that 

du/dy — ^ 03 as y — > 0, and Bj^ corresponds to vi^jK in the usual 

Incompressible case (cf. equation (47) and (W)), K being Von Karman' s 
Tiniversal constant. Approximately, within the range of moderate Mach 
nunibers, Bq^ might be taken to be independent of Mach niimber without 
serious error. 
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Some Experimental Results of Turbulent 
Boundary Layer 

% 

A brief examination of the velocity distribution in the incompres- 
sible case may be made to check with the logarithmic law predicted by 
similarity theory. A quite common form for turbulent-boundary-layer flow 
has been the l/7-power law. Recently, there are the careful measurements 
in the Rational Bureau of Standards^ and some others such as those by 
Hama (reference 10 ). Typical examples from these works at varioios 
Reynolds numbers have been plotted on semilogarithmic paper as fig- 
ure 1. It is found that outside the "viscous" layer near the wall, 
all the curves tend to show a more or less straight portion, of approxi- 
mately the same slope, the extent of the straight portion decreasing 
as the Reynolds number decreases.^ Such a phenomenon is in accordance 
with the expectation bn the basis of the similarity theory, because a 
smaller Reynolds number brings in the viscous effects to a greater degree, 
calls ing a merging of the higher- and lower-frequency parts of the energy 
spectrum. The slope was predicted to be proportional to u-j- and hence 
should change only slowly with Reynolds number, again verified by the 
experiments. For larger values of y/S, the velocity curve deviates from 
the straight portion on the semilogarithmic plot and may be approximated 
by a power law. In this part of the boundary layer, the turbulence is 
nearly isotropic and therefore the mean flow is not expected to follow 
closely the logarithmic relation. The empirical l/7-power law is also 
included and shows good agreement on the whole. 

For the compressible case, the works of Wilson (reference 3) and 
Ladenburg sind Bershader (reference 9 ) have been quoted before. The 
velocity distributions given are also plotted in figure 1 for comparison. 
It should be noted that in deducing the velocity distribution from their 
measurements, both made use of the isoenergetlc relation. Wilson’s 
measurements were done with a pitot tube and registered actually the 
Mach number variation within the boundary layer, computed from the 
stagnation pressures. A discussion was made on the error Introduced 
by the assumed isoenergetlc relation, with the conclusion that the 
maximum error amounts to about 3 percent for = 2.0. Row Ladenburg 
and Bershader measiired by means of an interferometer and, therefore. 


^As yet unpublished. The authors are obliged to Dr. Schubauer, 
chief of the Aerodynamics Section in RBS, for furnishing the data. 

^Hama's resvilts show more wavlness than those of RBS, partially 
because his Reynolds number is quite low, and are omitted in figure 1. 
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have recorded actually the temperature variation. By coincidence the 
two sets of experiments have some of the tests made under almost 
Identical, conditions, like the following: 

* 

Reference 9: = l43° C, = 2.3, Rg ~ 10^ 

Reference 8:^ = 150° C, Mj^ = 2.03, Rg - lO^ 

After correcting for the Mach number difference, one may estimate the 
temperature distribution in Wilson’s case from Ladenburg anfl Bershader's 
data and evaluate the velocity from the measured local Mach number. In 
correcting for the Mach n\miber difference, assume a modified isoenergetlc 
law to hold. 


= 1 + - 

V 

a being a constant, so that 
Then, letting subscripts refer to the respective experiments. 




and B 


- 1 


(%^)„ - (M 


L and B 


W) 


L and B 



^Table I, station 6, of reference 8. 
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The following resxilts are then obtained: 


y /5 

0.104 

0.217 

0.321 

0.477 

0.686 


.665 

.751 

.820 

.903 

.982 


1.31 

1.24 

1.20 

1.12 

1.05 

(u(ui)w 

.762 

.836 

.899 

.958 

1.003 


.769 

CO 

on 

CO 

• 

! 

. 

00 

00 

00 

.943 

.990 


where the subscripts Wff refer to the value given by Wilson iising the 
isoenergetlc law. It is Interesting to see that, at least for this 
particular case, the isoenergetlc law does give a very close approxi- 
mation as checked by the independent temperature measurement. However, 
whether the same agreement would result in other cases cannot be readily 
ascertained. 


Variation of Skin Friction with Eeyaolds 
and Mach Humbers 


In the present theory the skin friction was introduced as an 
external condition such as would influence the turbulence pattern 
through its effects on the mean distributions. One may stretch the 
theory a little to yield an eqxiation of the skin friction involving the 
choitLwi.se Reynolds number and the Mach number, as Prandtl and Von Kdrm&n 
had done for the incompressible case (references 11 and 12). A momentum 
thickness d for the entire boundary layer may be derived, which is 
related to the skin friction by 


T 


Pi ^1 


2 ^ 
dx 


The viscosity is introduced to form a length scale with the frictional 
velocity Uy and a chordwise Reynolds number appears. When the univer- 
sal consteints in the present theory are known, the skin- friction relation 
thus derived will also give explicit Mach nuniber effects. A point to be 
noted is that the velocity distribution should ref er to the frictional 
velocity u-j-, instead of u^^ at the outer edge y = 6. For simplicity 

it might be necessary to make approximations but there seems to be no 
great difficulty. 


Massachusetts Institute of Technology 

Cambridge, Mass., December 27, 1950- 



APPENDIX A 


SYMBOIS 




constants j defined by equations (33) and (39)^ 
respectively 


coefficients in expansions for cpj^ and cp£ 
defined in equations (5l) 


CfjCq 


'P^"v 


* * * 


• « « 


e ^ 
fl,f2 
SqjS'].^ * 

SqjS_2^^ ' 
k 

ki^kg 

^0 

m 


coefficients of skin friction and beat transfer, 
respectively, defined in equations (l8) 

specific beats at constant pressure and constant 
volume, respectively 

coefficients in expansions for Yu,l and '^■a,2> 
defined by equations (5^) and (55)^ respectively 


dilatation 


^ 8u ^ 8v ^ 

^8x 8y hz 


) 


correction functions for Macb ntmiber effect on 
logarithmic law, defined by equations (57) 

functions in expansion for G, defined by 
equations (4l) 

functions in expansion for i/G, defined by 
equations (te) 

constants defined by equations (56) 


coefficient of beat conductivity of fluid 
constants defined by equations (30) 

similarity scale of length 

rate of mass transfer at wall 
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P^Pl 


q 

-X- 

q 

Si,s2 

t 

u,v,w 




u* 




u 


mean pressure and mean pressure in free stream, 
respectively 

rate of heat tranter at irall 

rate of amount of heat generated in laminar sublayer 

constants defined by equations (30) 

time 

similarity scale of time for fluctuations 

velocity con^ranents in x-, y-, and z-directions, 
respectively 

similarity scale of velocity fluctuations 

modified mean velocity, defined by equation (16) 

magnitude of mean velocity at edge of lam i nar 
sublayer 

nondimensionalized modified mean velocity ^u = u^ju^) 


^1 


value of tl evaluated at outer edge of turbxaent 
boimdary layer • 


Uj, 

Ut 


value of u evaluated at a reference station within 
boundary layer 

frictional velocity 




x,y,z 


yr 

-A-i,-A2 

Bi,B2 


C 

P 

G 

Ml 


Cartesian coordinates; x, axis in direction of plate 
nnfl free stream; y, axis normal to plate; and z, 
axis parallel to leading edge of plate 

reference station in boimdary layer 

constants defined by equations (20) 

constants defined by equations (^5) 

constant defined by equations (3l) 
function defined by equation (2?) 

function defined by equation (35) 

Mach number of free stream 
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Q 

R* 

R5 

T 

U 

X 

^e,o’^e,v * • * 




3 


Pi 

7 

5 

6* 


e 

0 


®1^®2 




0 


( 0 )^"( 1 )^ 
(2)^ • • • 


constant defined by equations (20) 

gas constant in equation of state 

Reynolds number -based on botmdary- layer thickness 

temperature 

mean temperature at vaJ.1 

mean velocity component in x-direction 

function defining similarity scale of length, as 
in equation ( 32 ) 

distribution functions for velocity, defined 

by equations (-43) 

distribution functions for mean velocity, defined 
by equations ( 60 ) 

constants representing correlations in averaged 
equation of energy 



constant defined by equations ( 64) 

ratio of specific heats 
thickness of bovindary layer 
thickness of laminar sublayer 
rate of dissipation 

nondlmensionalized mean temperature, defined by 
equations ( 18 ) 

similarity scales of teng)erattire fluctuation 

functions in expansion for 0, defined by 
equations ^24) or ( 25 ) 
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1 ^ 

V 

p 

T 


To ^ Toz 


• • • 


momentum thickness of boundary layer 
coefficient of viscosity 
coefficient of kinematic viscosity 
density of fluid 
shearing stress 

shearing stress at the wall in direction of main 
stream and parallel to leading edge^ respectively 

functions in expansion for F, defined by 
equations (to) 


The subscript 1 denotes quantities in the free stream, Beirred 
quantities always represent mean values; primed quantities represent 
f luc tuat ions . 
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APPENDIX B 

DISSIPATION IN THE VISCOUS SUBLAYER 


Imagine the houndary layer to he idealized to consist of two parts; 
The major part is entirely txrrhulent, extending from the outer edge to 
a certain distance close to the wall. From there on the viscous effects 
take over completely, and a viscous sublayer is formed. In the viscous 
sublayer, the flow is laminar and free of turbulence. 

Let the thickness of the viscous s\iblayer be 8*. The energy 
equation within the sublayer then degenerates into 



(Bl) 


since all the correction terms vanish for the laminar flow with 
u = u(y), V = 0. The first term is the viscous dissipation; the 

dAX 

second is the heat conduction. Noting that u — =s = Constant, one 

dy 

gets by integration between 0 and 8*, 




T„u* - k 


dy 


y=0 


= 


(B2) 


assuming conduction to be conrparatively negligible at y = 8*. Obvi- 
ously TqU* Is the heat generated in the viscous sublayer by dlssipa- 
dT 

tlon, and -k — is the heat conducted into the fluid from the 

dy y=o 

wall. In the notation of equation (9)^ they are the qxiantlties q* 
and q, respectively. Therefore, q^ is the resxiltant heat exchange 

between the vlscovis and the turbulent layers. 

It is in 5 )ortant to estimate, at least approximately, the magni- 
tude u* for the present theory (cf. equation (l6)). To do so the 
first step is to estimate the thickness 8*. Usually the thickness 
of the sublayer is expressed by 
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<2.6*Ut 

V 


= h. 


a constant 


(B3) 


In the acttial case there is a transition layer where the viscous and 
turtmlent shears are of the same order of magnitude. The value of h 
in the incompressihle case varies from approximately 10 to 30, 
depending on whether the sublayer is taken to he strictly laminar and 
has, consequently, a linear velocity profile, or the sublayer is to 
consist of the entire thickness where noticeable departure from the 
turbulent velocity profile occurs. In the idealization adopted in this 
section, the transition layer is omitted and the boundary of the 
"viscoxis layer” should therefore lie somewhere within the actual tran- 
sition layer. By such reasoning, it seems that an average value of h 
of the order of 20 is probably adequate for the purpose. 

It may be seen that, with such a concept, the value of h ought 
to vary with the Mach ntmiber. For, one may regard the sublayer as owing 
its existence to the fact that the viscous shear becomes an appreciable 
part of the total shear. Suppose, then, that h is defined to satisfy 


du 

dy 


y=8* 


= T a 
0 


(b4) 


a being a constant less than unity. The quantity du/dy may be 
estimated by using the turbulent velocity distribution. For Instance, 
in the subsonic case if one tises equation (57), 


u. = — loge y + B + 
IC 




du 

dy 





Hence, by substituting into eqviation (b4). 
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■where u-^ = — u* (cf. equation (l6))- 
Incompresslhle case, there follows 


If ^ 


is the value for the 


^ = 1 + Mi^fi’(u*) + . . . (B5) 

hj^ - 


Since the right-hand side involves u*, successive approximations might 
he necessary for obtaining the solution. 

Having determined 6*, one may find u* by again applying the 
turbulent distribution. 
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